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LONG-TERM  GOALS 

This  programmatic  effort  has  three  long-tenn  goals.  The  first  is  to  simulate  bioacoustic  interactions 
within  and  near  individual  fish.  We  developed  a  methodology  that  combines  x-ray  CT  scans  with 
tissue  elasticity  measurements  and  finite-element  modeling  software  (WHACR).  This  technique  has 
provided  significant  insights  and  discoveries  regarding  toothed  whale  bioacoustics  (Cranford  et  ah, 
2008b;  Cranford  et  ah,  2008c),  and  now,  within  a  fish’s  head.  The  second  long-term  goal  is  to 
improve  and  refine  our  ability  to  measure  tissue  elasticity  in  samples  by  building  a  portable  device  to 
measure  physical  properties  from  tissue  samples.  The  third  and  final  goal  is  to  validated  the  finite 
element  models  and  add  more  functionality  to  the  WHACR  software. 

OBJECTIVES 

The  primary  objective  is  to  examine  simulations  of  the  otolith  organs  to  elucidate  their  motion 
patterns  in  response  to  acoustic  stimuli  from  different  directions  and  frequencies.  The  secondary 
objective  is  to  standardize  and  improve  our  ability  to  measure  tissue  properties;  primarily  Bulk 
Modulus  and  sound  speed  in  tissue  samples.  The  tertiary  objective  is  to  refine  the  vibroacoustic 
toolkit  (WHACR)  and  validate  these  finite  element  models. 
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APPROACH 


In  teleost  fishes  there  are  three  pairs  of  otoliths,  dense  masses  or  “stones”  of  calcium  carbonate,  each 
sitting  upon  a  patch  of  hair  cells,  all  contained  within  fluid  filled  sacs.  We  have  generated  high- 
resolution  CT-scans  of  the  White  Seabass  (Atractoscion  nobilis )  and  segmented  the  otolith  organs  from 
the  scanned  volumes  for  image  analysis. 

Otolith  organs  may  allow  fish  to  analyze  sound  frequency  and  direction,  but  the  mechanism  remains 
elusive.  It  is  necessary  to  account  for  the  ability  of  fish  to  discriminate  sounds  of  differing  frequency 
(reviewed  by  Enger)  (1973).  Sand  has  suggested  that  the  movement  patterns  of  the  otoliths  may  be 
frequency  dependent,  and  that  the  parts  of  the  macula  which  are  stimulated  may  depend  upon 
frequency  (Sand,  1974).  Further  studies  are  necessary  before  this  notion  can  be  confirmed. 

Fish  can  determination  the  direction  of  a  sound  source  (see  review  by  Sand  and  colleagues)  (Sand  and 
Bleckmann,  2007;  Sand  and  Karlsen,  1986).  One  obstacle  to  understanding  the  mechanisms  for  these 
abilities  is  our  lack  of  knowledge  about  the  movements  of  the  otoliths  themselves.  For  example,  do 
otoliths  show  simple  translation  back  and  forth  along  the  axis  of  sound  wave  propagation  or  are  the 
motions  more  complex?  Horner  (1980)  suggested  that  the  saccular  otolith  can  only  rock  within  the 
fluid-filled  sacs  and  upon  the  macula  because  of  constraints  from  proximal  structures  and  fibrous 
attachments  within  the  ear. 

In  order  to  determine  the  effect  that  otolith  shape  might  have  on  otolith  motion,  we  ran  two  sets  of 
simulations  using  WHACR,  our  numerical  analysis  software.  In  the  first  simulations  (Extracted 
Otoliths  Simulations),  we  extracted  the  shapes  of  all  six  otoliths  from  high-resolution  microCT  scans 
of  the  White  Seabass.  Those  “otoliths”  were  assigned  uniform  calcareous  material  properties, 
immersed  in  a  simulated  shear-soft  jelly,  and  exposed  to  different  stimulus  frequencies  and  directions. 

In  the  second  set  of  simulations  (Simplified  Otolith  Simulations),  we  compared  the  responses  of 
simplified  “otolith”  shapes  with  the  results  of  the  Extracted  Otoliths  Simulations.  In  this  second 
simulation,  we  used  two  simple  shapes,  a  “spherical  otolith”  and  a  “hemispherical  otolith.” 

One  crucial  step  in  building  a  finite  element  model  is  a  process  called  validation.  This  process  tests 
the  veracity  of  the  model  by  comparing  virtual  simulations  to  actual  experimental  results.  We  are 
currently  generating  experimental  results  that  will  be  compared  to  the  virtual  model  simulations. 

Our  first  attempts  to  develop  a  device  to  measure  tissue  elasticity  were  unsuccessful.  We  are  currently 
developing  an  alternative  methodology,  modifying  a  device  known  as  the  Linear  Skin  Rheometer 
(LSR)  as  the  central  component  (Hess  et  ah,  2006).  We  will  also  combine  this  device  with  a  means  to 
measure  sound  velocity  in  the  tissue  samples. 

WORK  COMPLETED 

We  have  examined  the  two  elementary  models  of  “otoliths”  (Extracted  Otoliths  Simulations  & 
Simplified  Otolith  Simulations)  and  report  that  the  otoliths  rock  in  response  to  planar  harmonic  waves 
from  different  directions  and  stimulus  frequencies.  We  have  drafted  two  manuscripts  that  report  on 
each  set  of  simulations.  These  reports  are  now  in  the  review  and  revision  process. 
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Development  of  new  tools  for  the  vibroacoustic  toolkit  (WHACR)  continues.  The  newest  tool  will 
provide  the  capacity  to  change  the  scale  of  element  size  during  simulations  within  the  model.  This  will 
allow  us  to  reduce  computational  overhead  for  large,  whole-animal  data  sets. 

A  procedure  for  extrapolating  the  scattered  or  total  harmonic  pressure  in  the  near  field  using  the 
Helmholtz  integral  was  implemented  in  the  vibroacoustic  toolkit.  This  tool  will  allow  us  to  calculate  a 
reverse  projection  of  the  head-related  transfer  function  (HRTF)  in  a  fraction  of  the  time  required  for 
the  forward  projection,  saving  time  and  computational  overhead.  We  have  tested  this  algorithm  on  the 
CT  dataset  from  a  live  bottlenose  dolphin  (acquired  from  the  US  Navy  Marine  Mammal  Program). 

We  will  also  calculate  the  forward  projection  of  the  sound  field  for  use  in  our  model  validation  efforts. 
The  calculated  HRTF  sound  field  will  be  compared  to  the  psychoacoustic  HRTF  measured  from  live 
bottlenose  dolphins.  These  measurements  with  live  dolphins  are  currently  underway  in  Paul 
Nachtigall’s  laboratory  on  Cocoanut  Island,  Hawaii. 

Challenges  and  Solutions 

Fish:  Anatomically,  it  is  important  to  know  the  precise  position  of  the  sensory  maculae  with  respect  to 
the  otoliths.  Discovering  this  anatomic  relationship  has  been  problematic  because  the  otoliths  are  so 
hard  that  they  cannot  be  effectively  sectioned  in  their  anatomic  context.  We  have  now  identified  a 
staining  method  that  will  allow  us  to  enhance  the  soft  tissues  in  the  CT  scans.  This  method  looks  very 
promising. 

Validation:  In  order  to  complete  the  model  validation  process,  we  need  to  acquire  a  fresh  bottlenose 
dolphin  specimen  from  the  stranding  network  at  the  University  of  North  Carolina  at  Wilmington.  They 
have  agreed  to  send  a  fresh  specimen  when  one  becomes  available.  Once  we  have  the  specimen  we 
can  quickly  move  forward  to  scan  the  specimen  and  run  the  requisite  simulations. 

Elasticity:  We  have  had  to  restart  our  efforts  to  build  a  device  that  can  measure  tissue  elasticity  from 
necropsy  samples.  This  was  a  setback  but  we  are  currently  working  on  a  new  trajectory  to  accomplish 
this  part  of  the  project.  The  LSR  device  will  be  built  and  tested  by  a  former  naval  officer  and 
advanced  engineering  graduate  student  as  part  of  completing  his  dissertation  research. 

Computer:  Running  vibroacoustic  simulations  using  finite  element  models  that  contain  millions  of 
elements  is  a  computationally  intensive  undertaking.  At  the  beginning  of  this  project,  we  purchased  an 
advanced  computer  that  contains  48  CPU  cores  and  128  GB  of  memory  to  satisfy  our  computational 
needs.  Recently,  this  computer  experienced  a  problem  and  had  to  be  sent  back  to  the  factory  to  be 
repaired.  We  eagerly  await  the  return  of  this  valuable  piece  of  equipment  so  that  we  can  continue  our 
modeling  work. 

RESULTS 

Extracted  Otoliths  Simulations 

We  collected  high-resolution  anatomic  data  from  a  few  small  (-21 -cm  total  length)  dead  Atractoscion 
nobilis  (Sciaenidae)  from  southern  California  by  means  of  a  micro-CT  scanner.  The  scan  data  was 
used  to  extract  high-fidelity  representations  of  the  otoliths  and  build  a  finite  element  model  (FEM)  by 
the  methods  and  tools  developed  by  Krysl  et  al.  (Krysl  et  al.,  2008)  and  Cranford  et  al.  (Cranford  et  ah, 
2008a;  Cranford  et  al.,  2008c).  This  FEM  model  allowed  us  to  investigate  the  dynamic  response  of 
fish  otoliths  to  incident  planar  acoustic  waves  whose  wavelengths  are  much  longer  than  the  dimensions 
of  the  otoliths.  The  otoliths  are  modeled  as  embedded  in  a  shear-soft  fluid-like  jelly.  The  simplified 
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model  does  not  currently  include  any  other  structures,  such  as  the  nearby  cranial  bones  nor  influences 
of  the  swimbladder.  The  model  space  was  simulated  with  two  different  sinusoidal  signals  (200  and 
400  Hz)  from  several  different  directions  with  respect  to  the  fish. 

Shear  forces  result  from  the  relative  motion  between  the  otolith  surfaces  and  the  shear-soft  jelly  that 
surrounds  them  in  the  model  space.  The  results  show  that  the  400-Hz  simulation  produces  greater 
shear  values  (due  to  larger  displacements  of  the  otoliths),  particularly  in  the  dorsoventral  dimension, 
than  does  the  200-Hz  signal  of  the  same  magnitude  and  direction.  Similarly,  changing  the  direction  of 
the  acoustic  stimulus  also  produces  altered  patterns  of  shear  forces  acting  upon  the  surfaces  of  the 
otoliths.  These  FEM  simulations  produced  intriguing  results,  suggesting  that  frequency  and  direction 
may  be  encoded  by  the  unique  rocking  motions  of  the  otoliths.  If  this  rocking  behavior  holds  true  for 
actual  otoliths  then  we  will  have  discovered  a  basic  mechanism  for  hearing  acuity  in  fish. 

Simplified  Otolith  Simulations 

The  intriguing  success  of  the  Extracted  Otolith  Simulations  catalyzed  an  investigation  into  whether  the 
simplest  “otolith”  shapes  would  produce  similar  intriguing  results.  In  essence,  do  “spherical”  or 
“hemispherical”  shapes  (scatterers)  produce  angular  oscillations  due  to  torque  in  the  presence  of 
relatively  long  wavelength  acoustic  stimuli?  These  simple  scatterers  can  be  considered  as 
“abstractions”  of  otoliths  because  they  are  without  the  distinctive  boat-shape  and  sculpting  found  on 
actual  otoliths,  features  which  are  often  used  as  keys  to  species  identification.  In  these  Simplified 
Otoliths  Simulations  we  considered  a  progressive  planar  harmonic  wave  in  an  acoustic  fluid  with 
selected  mass  density,  speed  of  sound,  and  frequency 

The  acoustic  waves  impinge  upon  a  stiff  homogeneous  scatterer  (the  simplified  “otolith”)  of  arbitrary 
shape,  whose  characteristic  dimensions  are  all  much  smaller  than  the  wavelength  of  the  incident 
acoustic  wave.  In  the  case  of  a  “spherical  otolith,”  uniform  pressure  across  the  symmetrical  sphere  in 
the  acoustic  fluid  does  not  generate  an  accelerating  torque  on  the  scatterer.  However,  the  case  of  the 
hemispherical  scatterer  is  different.  There  is  a  dynamic  torque  experienced  by  the  hemispherical 
scatterer.  The  numerical  analysis  indicates  that  the  X  and  Z  components  of  the  torque  are  identically 
zero  (these  are  the  directions  of  the  nonnal  to  the  propagating  sound-waves,  and  the  direction  of  the 
axis  of  the  symmetry  of  the  scatterer).  However,  the  hemispherical  scatterer  oscillates  (rotationally) 
about  the  Y axis  (in  this  direction  is  parallel  to  the  propagating  wave  fronts).  It  rocks! 

Furthennore,  for  other  non-spherical,  asymmetrical,  shapes,  the  scattered  pressure  will  generate  a 
dynamic  torque  which  will  result  in  angular  motion  (rocking)  in  accordance  with  the  specific 
parameters  and  initial  conditions.  This  principle  suggests  that  all  such  “asymmetric”  shapes,  like 
otoliths  will  experience  rocking  in  the  presence  of  relatively  long  wavelength  progressive  wave 
stimuli. 

We  have  drafted  and  submitted  a  manuscript  about  these  simplified  synthetic  scatterer  results  and  it  is 
currently  in  the  process  of  revision,  pending  changes  for  calculations  of  flow  characteristics. 

Anatomy  of  Otolith  Organs  and  Sensory  Maculae 

The  anatomic  relationships  between  the  otoliths  and  the  sensory  maculae  are  the  basis  for  the  input  to 
the  fish  central  nervous  system.  Discovering  these  relationships  has  been  an  intractable  problem 
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caused  by  the  hardness  of  the  otoliths  near  the  delicate  nervous  tissue.  Fortunately,  we  have  found  a 
pre-scan  Iodine  staining  technique  developed  by  Metscher  and  colleagues  (Metscher,  2009a;  b)  that 
will  make  these  relationships  visible  after  CT  scanning. 

In  addition,  we  have  salvaged  the  sensory  maculae  from  five  White  Seabass  specimens  so  that  we  can 
obtain  a  map  of  the  hair  cell  directions.  The  hair  cell  maps  are  currently  being  produced  in  the 
laboratory  of  Dr.  Art  Popper  in  Maryland.  The  line  scale  motions  of  the  otoliths  combined  with  the 
hair  cell  maps  may  provide  clues  to  the  means  by  which  fish  can  discriminate  the  frequency  and 
direction  of  a  sound  source. 

Model  Validation  and  Head-related  Transfer  Function 

The  experimental  design  for  the  HRTF  experiments  with  live  dolphins  in  Hawaii  has  been  finalized 
and  data  collection  has  begun.  That  portion  of  the  study  should  be  complete  by  the  end  of  the  year. 

The  construction  of  the  model  that  we  will  use  for  comparative  purposes  has  started  with  the  CT 
dataset  from  the  live  bottlenose  dolphin  acquired  from  the  US  Navy  Marine  Mammal  Program.  We 
will  use  this  dataset  until  we  obtain  a  fresh  postmortem  specimen  from  UNCW  and  can  generate  a 
higher  resolution  model. 

IMPACT/APPLICATION 

The  Extracted  Otoliths  Simulations  and  the  Simplified  Otolith  Simulations  produced  similar  “rocking” 
motions  from  unsymmetrical  “otoliths.”  This  consistency  hints  at  a  potential  principle  of  otolith 
motion  that  might  be  common  to  all  teleost  fish.  Future  iterations  will  include  additional  anatomic 
components  within  the  model. 

RELATED  PROJECTS 

This  project  is  an  outgrowth  of  the  methodology  we  have  developed  over  the  past  eight  years 
(Cranford  et  ah,  2010;  Krysl  et  al.,  2008;  Krysl  et  al.,  2006;  Soldevilla  et  al.,  2005).  We  are  currently 
using  the  same  basic  methods  to  study  interaction  between  toothed  whale  anatomy  and  selected 
sounds. 

PUBLICATIONS 
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software.  Report  to  University  of  California.  San  Diego,  p  9. 

Krysl  P,  Trijoulet  V,  Cranford  TW.  (In  press).  Validation  of  a  vibroacoustic  finite-element  model  using 
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Life.  New  York:  Springer  Science+Business  Media,  LLC. 
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Mammals. 
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